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Compactness of Cold Gas in High-Redshift Galaxies 
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ABSTRACT 

Galaxies in the early Universe were more compact and contained more molecular gas than 
today. In this paper, we revisit the relation between these empirical findings, and we quantita- 
tively predict the cosmic evolution of the surface densities of atomic (HI) and molecular (H2) 
hydrogen in regular galaxies. Our method uses a pressure-based model for the H 2 /HI-ratio of 
the Interstellar Medium, applied to ~ 3 • 10 7 virtual galaxies in the Millennium Simulation. We 
predict that, on average, the Hi-surface density of these galaxies saturates at Shi < 10 M© pc~ 2 
at all redshifts (z), while H2-surface densities evolve dramatically as Sh 2 00 (1 +£) 2 ' 4 - This scal- 
ing is dominated by a oc (1 + z) 2 surface brightness scaling originating from the oc (1 + z)~ l 
size scaling of galaxies at high z- Current measurements of Xh 2 at high z, derived from CO- 
observations, tend to have even higher values, which can be quantitatively explained by a 
selection bias towards merging systems. However, despite the consistency between our high-z 
predictions and the sparse empirical data, we emphasize that the empirical data potentially 
suffer from serious selection biases and that the semi-analytic models remain in many re- 
gards uncertain. As a case study, we investigate the cosmic evolution of simulated galaxies, 
which resemble the Milky Way at z = 0. We explicitly predict their HI- and H2 -distribution at 
Z= 1.5, corresponding to the CO-detected galaxy BzK-21000, and at z = 3, corresponding to 
the primary science goal of the Atacama Large Millimeter/submillimeter Array (ALMA). 

Key words: ISM: atoms - ISM: molecules - ISM: clouds - galaxies: evolution - galaxies: 
high-redshift. 



1 INTRODUCTION 

Gala xies were more compact in the ea r ly Un i verse than to- 
day jBouwens et all 12004 iTruiillo et all 120061 : iBuitrago et aD 
2008). This empirical feature is probably driven by an increase 
in the volume-to-m as s ratio of dark haloes with cosmic time 
dGunn & Gottlll972h . which dictates the d ensity evolution of galax - 
ies by the transfer of ang ular momentum (Fall & Efstathiou 1 9801) . 
We have recently argued ( Obreschkow & Rawlingsll2009ah that the 
density evolution of galactic discs implies a systematic pressure 
change, which causes a dramatic decline in the mass ratio between 
molecular (H 2 ) and atomic (HI) hydrogen with cosmic time. Us- 
ing the Millennium Simulation, we showed that this decline in the 
H 2 /HI-ratio simultaneously explains (i) the observation s of large 
molec ular masses in ordinary galaxies at z = 1.5 dPaddi et aD 
2008), (ii) the weak cosmic evolution of the Hi-density £2 H i in - 
f erred from damped Lyman-of systems (Prochaska et al ] 120051) . 
and (iii) the history of star formation inferred from ultravio- 
let, far-infrared, submillim eter, and radio continuum observations 
(iHopkins & BeacomEooih . 

In this paper, we quantitatively predict the cosmic evolution 
of the surface densities of HI and H 2 in a large sample of regular 
galaxies. Section [2 explains our physical model and the numerical 
simulation. Section [3] presents our predictions for the average cos- 
mic evolution of the HI- and H 2 -distributions in a broad sample of 
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Figure 1. Surface densities of galactic cold gas as a function of galactic 
radius. Lines show the predictions of Eqs. {T][2} f° r ^hi(^) (solid), I<u 2 (r) 
(dashed), and Zin(r) + ^n 2 (r) (dotted) for different values of R^ oV Dots 
represent the observed counterparts iLerov et al. 2008) for E ffl (r) (filled) 
and Sn 2 (r) (empty), inferred from CO(2-l) observations. 



galaxies, while Sectional specifically focusses on the evolution of 
HI- and H 2 -distributions in Milky Way (MW)-type galaxies. Sec- 
tion [5] gives a brief summary and outlook. 
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2 SIMULATION OF E ffl AND E H2 IN GALAXIES 



Several computational studie s (e.g. [Co mbes 



et al.1 

ivtsovl 



1999 



2008; 



Greve & Sommer-Larsenl 120081 : iRobertson & Kravtsa.. _ 

Dobbsll2008h lave investigated the distribution of molecular g; 
in regular galaxies and revealed pivotal connections between 
the distribution of mo l ecules and s tar formation. I n parti cular, 
IRobertson & Kravtsovl d2008h and iKrumholz etall 120091 (see 
also Elmegreer] Il993 ) demonstrated that the empirical relation 
betw een the local H 2 /H I -ratios and cold gas pre s sure, measured 
by iBlitz & Rosolowskvl d2006h and iLerov et al.l d2008h . can be 
approximately reproduced by SPH/Af-body simulations, which 
include a model for interstellar radiation. 

While those simulations are effective tools for the study of the 
cold gas in individual galaxies, they are computationally too ex- 
pensive to be included in "cosmological simulations", i.e. simula- 
tions of representative samples of galaxies with resolved merger 
histories in an expanding model-universe. However, cosmologi- 
cal simulations of HI- and H 2 -distributions are required for the 
design and analysis of high-z cold gas surveys possible with fu- 
ture telescopes, such as the Square Kilometer Array (SKA), the 
Large Millimeter Telescope (LMT), and the Atacama Large Mil- 
limeter/submillimeter Array (ALMA). To circumvent the current 
computational bottleneck, we decided dObreschkow et alj|2009h to 
adopt a semi-analytic simulation of millions of galaxies and to es- 
timate the HI- and H 2 -distributions of these galaxies using an an- 
alytic model. We shall first explain this model and then the semi- 
analytic simulation. 

Our model for the distributions of HI and H 2 in regular galax- 
ies repre sents a rami f ication of the H 2 /HI-pressure relation mea- 
sured by ILerov et all d2008h , combined with two assumptions for 
the distribution of cold hydrogen. The first assumption is that all 
galaxies carry their cold gas in flat discs, such as i s observed in vir- 
tually all regul ar galaxies (s ee lLerov et al.l d2008h for spiral galax- 



ies at low z. Young j2002h for elliptical galaxies at low z, and 
alJdZI 



e.g. i Tacconi et al. (2006) for galaxies at high z). The second 
sumption is that the surface density of cold hydrogen, E ffl (r) + 
E H o(r), follows an exponential profile in the radial coordinate r (e.g. 
ILerov ctalJl2008h. Based on these assumptions we have shown 
(Obreschkow et al. 2009) that the surface densities of HI and H 2 are 



2in(r) 
2 H2 (r) 



£ H exp(-r/r disc ) 
1 +^ ol exp(-1.6r/r disc ) 
t n R c mol exp(-2.6 r/r disc ) 
1 +/^exp(- 1.6 r/rdisc) 



(1) 

(2) 



where r disc is the exponential scale radius of the cold gas disc, 
S H = (Mm + Mn 2 ) I \2nr 2 dhc ) is a normalization factor, and R° mol = 
^Ho(Q)/^m(Q) is the cent ral H 2 /HI-ratio. R c mol can be approximated 
( Obreschk ow et alJl2009h from r disc , the disc stellar mass M^ 1SC , and 
the total (HI+H 2 +He) cold gas mass M gas as, 



*moi = [K r d "it c M gas (M gas +0.4 Mf «)]° 



(3) 



where K - 1 1.3 m 4 kg -2 is an empirical constant. 

If the radii and surface densities in Eqs. CJO are normalized 
to r disc and E H , then E H i(r) and E H2 (r) become unique functions of 
R c moV Fig. \T\ shows E ffl (r) and £ H2 (r) in these normalized coordi- 
nates for three different values of R° „ chosen to match those of 

mol' I 1 

nearby galaxies with measured HI- and CO-densities (ILerov et al.1 
2008): NGC 4214, a star-forming, Hi-rich, dwarf galaxy; NGC 
5055, a massive spiral galaxy with similar HI- and H 2 -masses; 
NGC 5194, a barred, H 2 -rich spiral galaxy. The scale lengths r disc 



and the densities E H of these galaxies were estimated by fitting 
exponential functions to £ HI (r) + ^n 2 ( r ) f° r eacn galaxy (NGC 
4214/5055/5194: R c mol = 0.1/11/34, r disc = 2.3/5.1/2.6 kpc, S H = 
10/56/139 M©pc -2 ). The good fit between the model and the ob- 
servations demonstrates the validity of the pressure-based model 
for the H 2 /HI-ratio. 

Eqs. (0] O make a list of predictions, which can be 
seen as direct consequences of the H 2 /HI-pressure relation: (i) 
the H 2 /HI-mass ratio o f the whole galaxy increases with R° mol 
dObreschkow et al.l20 09): (ii) H 2 -rich galaxies, i.e. R c mol > 10, have 
annular Hi-distributions with a density drop towards the centre; 
(iii) if E ffl (r) and ^n 2 (r) cross, they do so close to the maximum 
of SHi(r). 

To evaluate E ffl (r) and ^n 2 (r) using Eqs. ([TJQ, we requir e 
an estimate of r disc , Mf lsc , and M gas . In lObreschkow et all |2009), 
we therefo re adopted the cosm ological galaxy simulation per- 
formed by iDe Lucia & Blaizotl d2007|) on the dark-matt er skele- 
ton of the Millennium Simulation ( Sprin gel et al ] [2005b . In this 
"semi-analytic" simulation, galaxies were represented by a list of 
global properties, such as position and total masses of gas, stars, 
and black holes. These properties were evolved using simplistic 
formulae for mechanisms, such as gas accretion by infall and merg- 
ers, star formation, gas heating by supernovae, and feedback from 
black holes. The free parameters in this model were tuned mostly t o 
optical observations in the local universe (see lCroton et al.1 120061 
such as the joint luminosity/colour/morphology distribution of op- 
tically observed low-redshift galaxies. However, no measurements 
of HI and H 2 were used to adjust the free parameters. 

The semi-analytic simulation resulted in a catalog listing the 
properties of ~ 3 • 10 7 model-galaxies at 64 cosmic time steps. 
In these galaxies, the cold gas was treated as a single compo- 
nent, hence masking the complexity of atomic and molecular 
phases. We therefore post-processed this simulation using Eqs. (0]- 
[3]) to evaluate HI- and H 2 -distributions for every galaxy. The 
resul ts of the e merging HI and H 2 simulations were presented 
in lObreschkow et al.1 d2009h . They well match the HI- and H 2 - 
mass functions, mass-diameter relations, and mass- velocity rela- 
tions observed in the local Universe. The high-z predictions are 
roughly consistent with the sp arse cold gas detections at z > 
dObreschkow & Rawlingsll2009ah . 

The limitations and uncertainties of this model for HI and 
H 2 at low z and high z were discussed in detail in Section 6 of 
Obreschkow et al . (2009). Some crucial points shall be emphasized 
here: (i) The predicted space density of HI, Q H i, undershoots the 
inferences of damped Lyman-a systems (DLAs) at z > 1 by about 
a facto r 2 (see Fig. 2 and Section 4 in lObreschkow & Rawlings! 
2009a). We discussed possible reasons for this discrepancy, but 
it should be stressed that an offset by a factor 2 is rather small , 
given that the semi-analytic models of IDe Lucia & Blai zot (2007) 
have never been tuned to Hi-data, but only to local optical data, 
(ii) One might worry whether the H 2 /HI-pressure relation, assumed 
constant in our model, coul d undergo a significant ev olution with 
redshift. In Section 6.3 of lObreschkow et al.l 120091), we argued 
that this effect is reasonably small because the empirical H 2 /HI- 
pressure relation seems to implicitly include changes in the photo- 
dissociating UV-radiation field, and th e scatter in H 2 /HI shows n o 
significant correlation with metallicity dBlitz & Rosolow skv 2006), 
although the mean metallicities of the galaxies in the sample dif- 
fer by a factor 5. This factor is on the order of the typical metal- 
licity evolution of an individual massive galaxy between redshift 
z = 10 and z = 0, and hence the cosmic evolution of metallicity 
seems to be a minor concern (for details and references see Sec- 
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Figure 2. Solid and dotted lines show the simulated cosmic evolution of 
the surface densities <S H ^ HM > and (X™ HM > with 1-cr scatter. These lines 
represent interpolations between 51 discrete redshifts between z = and 
z= 10. At each redshift, the average value relies on ~ 10 6 — 3 • 10 7 galaxies 
(the precise number depends on z). Dashed lines are the power law fits of 
Eqs. El 0. Points show the empirical data of Table [T] 



tion 6.3 in lObreschkow et aP l2009). (iii) Eq.[T] suggests that £ ffl (r) 
falls of exponentially at large galactic radii, while some observa- 
tions have suggested that it varies rather as 1/r, proportion al to the 
dark matter density derived from flat rotation curves (e.g. iBosmal 
1 19811 : iHoekstra et ail 2001). The detailed Hi-maps and CO-maps 
presented by iLerov et al.l d20Q8h for a sample of 18 nearby galax- 
ies nevertheless reveal that all these galaxies carry most of their 
cold gas mass inside a radius where the angularly averaged HI- 
distribution is consistent with an exponential profile. Therefore, we 
consider the exponential profile as a safe assumption for the pur- 
pose of splitting the cold gas mass between HI and H 2 . 

We use sample of 18 nearby galaxies, for which ILerov et al.1 
(2008) derived inclination-corrected column density maps of HI, 
H 2 , stars, and SFRs (not used here). 



3 COSMIC EVOLUTION OF 2 ffl AND Z H2 

In this section, we shall investigate the average cosmic evolution 
of the surface densities of HI and H 2 in the simulated galaxies 
(see Section[2]). We define the "average" (...) as the (M ffl + M H2 )- 
weighted geometric average over all galaxies with M m + M H2 > 
10 8 M o . The threshold M m + M H2 = 10 8 M o approximately marks 
the completeness limit of the simulation and the mass-weighting 
ensures that massive galaxies with low space densities contribute 
significantly to the average. 

For each galaxy we define the Full-Width-Half-Maximum 
(FWHM) surface densities, £™hm and e fw™ as me average sur _ 
face densities inside the radii r, where EmW and £ H2 (r) reach 
50% of their maximal value. The simulated cosmic evolution of 
<s fwh M) and ( £ fwhm ) is displayed in Fig. [2] as solid lines. Both 
<e fwhm )(z) and <s fwhm )(z) are weU fitted by power _i aws (dashed 



lines in Fig. [2]) 

<^r M /[M pc- 2 ]>fe) 

(E^ HM /[M pc- 2 ]>(z) 



4.7(1 +z) 02 , 
7.1 (1+z) 2 - 4 . 



(4) 
(5) 



The simulated Hi-density remains roughly constant reaching 
maximal values around 8-9 M© pc~ 2 , consistent with th e observed 
HI- saturation level in nearby galaxies dBigiel et al.l2008l) . This cos- 
mic evolution of is a pure prediction, with no observational 
counterpart at z > to-date. 



The simulated H 2 -densities can be compared to CO-data. The 
data points in Fig.[2](see also Table[T} represent 12 local and 1 1 dis- 
tant galaxies with observational estimates of E^ WHM derived from 
n 

resolved CO-mapo To our knowledge, none of the high-redshift 
galaxies are gravitationally lensed, however they certainly repre- 
sent a highly-biased sample (see below). 

At z = 0, the observed average of S R ^ HM , weighted by the 
cold gas masse s and the space densities drawn f rom the cold gas 
mass function dObreschkow & Rawlingsl l2009bh . is <£™ HM > = 
7 + 1 M pc" 2 , consistent with the simulated value of (S™ JM ) = 
7.9M pc" 2 . 

At z > 0, the empirical values of (E^™) show a clear trend 
to increase with redshift z. We must worry whether this is a pure 
selection effect, merely reflecting t he relation bet ween sensitivity 
and redshift (similar to Fig. 5 in lGreve et al.l l2005). To our knowl- 
edge, all the high-redshift CO- sources considered here have origi- 
nally been selected from single-mirror surveys (submillimeter, IR, 
optical), which are flux limited rather than surface brightness lim- 
ited. Unfortunately, it is not clear what fraction of all flux selected 
sources has been detected in the follow-up interferometric surveys, 
but it is likely that the present sample of spatially resolved high- 
redshift CO-sources is biased towards high fluxes (i.e. large H 2 - 
masses) rather than towards high surface brightness (i.e. high H 2 - 
surface densities). In general, more massive disc galaxies tend to 
have lower surface densities than less massive ones. Therefore, the 
selection bias of a flux limited sample tends to bias high-redshift 
disc galaxies towards lower H 2 - surface densities. If the observed 
systems were single discs, one would thus expect the real relation 
between H 2 - surface densities and redshift to be even steeper than 
suggested by the data points in Fig. [2 On the other hand, the CO- 
detected galaxies at z > seem to be heavily biase d towards sys- 
tems subjected to major mergers (iTacconi e t al. 2006). On average, 
major mergers decrease the specific angular momentum, hence de- 
creasing r disc and increasing S^hm 

We should therefore compare the CO-measurements at z > 
exclusively with the simulated galaxies with major mergers in their 
cosmic evolution. In the semi-analytic model, the evolution of the 
angular momentum during a merger is directly computed from the 
underlying Af-body Millennium Simulation. We can therefore ex- 
pect that, on average, mergers decrease the specific angular mo- 
mentum, and hence increase the cold gas densities towards the 
high obse rved densities sh own in Fig. [2] Within the semi-analytic 
model dCroton et al.l l2006\ where stellar bulges arise during merg- 
ers, galaxies resulting from major mergers are "elliptical" galaxies, 
defined as the obje cts with a bulge- to- tota l mass ratio larger than 
0.4 (see eq. 18 in bbreschkow et alJl2009h . The average (S^ HM > 
and the corresponding 1-cr scatter of the simulated elliptical galax- 
ies are represented by the green lines in Fig. [2l and provide a much 
better fit to the observational data. 

Despite the reasonable agreement between the simulated and 
the empirical data for (£™ HM ) and (E™ HM ), we emphasize that 
there are many subtle selection effects in the empirical data 
(e.g. flux limited sample versus surface brightness limited sam- 
ple, Eddington bias) and that the semi-analytic models remain in 
many regards simplistic and uncertain (e.g. no geometrical model 



1 The adopted CO-to-H 2 conversion is Xqo = 4 Kkms 1 pc 2 at z = 
(Table Q] top), co nsistent with observations in local quiescent galaxies 



LU top), co n; 

iLerov etal. 2008); and X C o = 1 Kkm s _1 pc 2 for sources at z > (Tabled 
bottom), consiste nt with ob servations in Ultra Luminous Infrared Galaxies 
iDownes & S olomon 199^). 
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Object 


Type 


z 


CO- 
line 


log(I™ 
[M pc- 2 ]) 


NGC0628 


Sc 


0.0 


2-1 


-0.4±0.1 (a) 


NGC3198 


SBc 


0.0 


2-1 


0.7±0.1 (a) 


NGC3184 


SBc 


0.0 


2-1 


0.6±0.1 (a) 


NGC4736 


Sab 


0.0 


2-1 


1.2±0.1 (a) 


NGC3351 


SBb 


0.0 


2-1 


0.1 ±0.1 (a) 


NGC6946 


SBc 


0.0 


2-1 


2.4±0.1 (a) 


NGC3627 


SBb 


0.0 


2-1 


1.6±0.1 (a) 


NGC5194 


SBc 


0.0 


2-1 


1.1 ±0.1 (a) 


NGC3521 


SBbc 


0.0 


2-1 


0.8±0.1 (a) 


NGC2841 


Sb 


0.0 


2-1 


-0.1 ±0.1 (a) 


NGC5055 


Sbc 


0.0 


2-1 


1.8±0.1 (a) 


NGC7331 


SAb 


0.0 


2-1 


0.8±0.1 (a) 


BzK-21000 


Galaxy 


1.5 


2-1 


2.1 ±0.2 (b) 


SMM J123549+6215 


SMG 


2.2 


3-2 


3.6±0.1 (c) 


SMM J163650+4057 


SMG 


2.4 


3-2 


3.3 ± 0.2 (c) 


SMM J163658+4105 


SMG 


2.5 


3-2 


3.5 ± 0.2 (c) 


SMM J123707+6214SW 


SMG 


2.5 


3-2 


2.4 ± 0.2 (c) 


SMM J123707+6214NE 


SMG 


2.5 


3-2 


2.6 ± 0.2 (c) 


BRI 1335 0417 


QSO 


4.4 


2-1 


2.9 ± 0.2< d > 


BRI 1202-0725 north 


QSO 


4.7 


2-1 


3.6 ± 0.2< d > 


BRI 1202-0725 south 


QSO 


4.7 


2-1 


4.0 ± 0.2 (d) 


Jl 148+5251 


QSO 


6.4 


3-2 


> 3.6< e) 



Table 1. Local (top) and distant (bottom) galaxies with spatially resolved 
CO-detections, and corresponding estimated H 2 -surface densit ies ^fwhm 
(a) From radial H2 -density profiles given in iLerov et al. ( 200 8|); (b) from 
CO-fluxes and velocit y peak separation m easured bv iDaddi et al J |2008); 
(c) from Table 1 in iTacconi et al.l (2006), for the submillimeter galaxy 
(SMG) SMM J 163650+4057 the radius was averaged between semi-major 
and semi-minor axis, for SMM J123549+6215 the CO-flu x was deter- 
mined from CO(3-2) and the radius from CO(6-5); (d) from Carilli et al. 
(2002), BRI 1335 0417 probably has two non-resolved sub-components; 
(e) IWalteretafl ( l2004l) . 



of mergers, mass resolution limitations, no ram-pressure stripping). 
Clearly, the advent of telescopes, such as the SKA, the LMT, and 
ALMA, will prove invaluable to refine these models. 

What are the reasons for the different cosmic evolutions of 
E ffl and E H2 ? According to Eqs. CHS, these evolutions can be 
understood from the average cosmic evolution of r disc , M gas , and 
M s . Most of the massive galaxies in the semi-analytic simulation 
are gas-dominated at z > 1 and their cold gas masses remain 
roughly constant with cosmic time due to a self-regulated equilib- 
rium between the net cold gas accretion and star formation. Thus, 
for any given galaxy, the cosmic evolution of the value R° mol and 
the functions 2 ffl (r) and ^n 2 (r) is essentially dictated by the evo- 
lution of the scale radius r disc . Assuming a similar specific angula r 
momentum for the galaxy and its halo dFall & Efstathiou| [l980), 
r disc is expected to evolve proportionally to the virial radius r vir 
of the halo. For a spherical halo of mass M vi r, the latter scales a s 
r 3 ir oc Mvir/[n m (l + z) 3 + Q A ] (flat universe. iGunn & Gotdfl972h . 
and hence for a fixed M vir , 



r d isc * [H m (l + z) 3 + Q A 



-1/3 



(6) 



where Q m (here 0.25) and Q A (here 0.75) denote the normalized 
space densities of matter and vacuum energy. 

At high z, Eq. J6]) reduces to r disc oc (1 + z)~ l , consistent wit h 
observations in the Hubble Ultra Deep Field teouwens et al]2 004). 
If the cosmic evolution of M gas and M s is neglected, then r disc oc 
(1 + zT l implies R c mol oc (1 + z ) 32 (see Eq.[3j and l u oc (1 + z ) 2 . 
Yet, for R c mol > 5/3, the maximum of E ffl (r) is given by E™ ax = 
0.516 E H ^ ol " 5/8 (see lObreschkow et alJExM) and hence E™ x oc 



(1 +z) 2 (1 +z) 2 = const. In other words, the surface density of HI is 
expected to show little evolution with redshift, consistent with the 
numerical fit for <S™ HM >(z) of Eq. {§. On the other hand, E H2 (r) 
in Eq. reduces to E H exp(-r/r disc ) V r > 0, if R c mol » 1 (i.e. z » 
1). Hence, if the cosmic evolution of M gas is negligible, £™ x and 
£™hm are predicted to scale as oc E H oc (1 + z) 2 . A comparison 
to Eq. © confirms that, within our model, the cosmic evolution of 
Eh 2 is largely explained by the size-evolution of galaxies, while the 
evolution of M gas plays a minor role, accounting for an additional 
oc (1 + z) 0A scaling. 



4 EVOLUTION SCENARIO FOR THE MILKY WAY 

We shall now investigate the cosmic evolution of the HI- and 
H 2 -distributions in MW-type galaxies. By definition, a simulated 
galaxy at z = is called a "MW-type" galaxy, if its morphological 
type is Sb-Sc and if it matches the stellar mass M s , the Hi-mass 
Mm, the H 2 -mass M H2 , the Hl-half-mass radius r H alf , and the H 2 - 
half-mass radius rJl alf of the MW within a factor 1.3. This factor 
roughly matches the empirical uncertainties of the MW data in Ta- 
ble O Within this definition, the simulation contains 2 • 10 3 MW- 
type galaxies at z = 0. At z > 0, we define MW-type galaxies as 
those objects, which are the most massive progenitors of a MW- 
type galaxy at z = 0. Most MW-type galaxies accrete about half 
of their mass in a time-interval corresponding to the redshift range 
z = 2 - 10. They typically undergo a series of minor mergers, al- 
lowing the build-up of small bulges. 

To accommodate the nature of the empirical data (e.g. poor 
data for gas at the MW centre), we here consider the half-mass 
radii r^ alf and rj^ alf , and the enclosed average surface densities E^ alf 
and lt lf , rather than FWHM values. 

In Table [2] we explicitly present the predicted average gas 
masses, gas radii, and gas densities of the 2 • 10 3 MW-type galax- 
ies at three specific redshifts: z = 0, corresponding to the MW it- 
self; z = 1.5, corresponding to the galaxy 5^-21000, which is the 
only "relatively ordinary" galaxy (i.e. a borderline ultra luminous 
infrared galaxy for a "normal" star formation efficiency) at higher z 
detected and spatially resolved in CO-emission to-date; z = 3, cor- 
responding to the primary science goal of ALMA, i.e. the detection 
of MW-type galaxies at z = 3 in less than 24 hours observation. 
Fig. [3] illustrates the average density profiles S ffl (r) and £h 2 (/) at 
these three redshifts. 

The cosmic evolution of the average HI- and H 2 -properties of 
the simulated MW-type galaxies can be characterized as follows: 
(i) the mass ratio M H2 /M H i approximately varies as (1 + z) L6 , iden- 
tical to the space density evolution £2 H2 /^hi oc (l+z) 1 - 6 predicted in 
lObreschkow & Raw lings (2009a); (ii) the surface density of HI re- 
mains approximately constant, while the density of H 2 increases 
by a factor 30 from z = to z = 3; (iii) both the Hi-radius and the 
H 2 -radius shrink by a factor 4-5 from z = to z = 3; (iv) the 
Hi-distribution is more annular at high z (see Fig. [3]), resembling 
the central Hi-deficiency seen in H 2 -rich local g alaxies. 

The comparison of the galaxy £z^-21000 (IDaddi et al.ll2008h 
and the MW-type galaxies at z = 1.5 (see Table [2]) suggests that 
the H 2 -density of BzK-21000 is consistent with that of a typical 
MW progenitor, while the total H 2 -mass of BzK-2\000 may be a 
few times higher. Considering that the MW lies at the lower end 
of intermediate mass spira l galaxies (stell ar masses M s in the range 
3 • 10 10 - 3 • 10 11 M^: lFlvnn et"aD l2006h . the H 2 -distribution of 
Z?z^-21000 appears to be typical for the progenitors of intermediate 
mass disc galaxies. 
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Figure 3. Cosmic evolution of the average profiles Sin(r) and £H 2 ( r ) cal- 
culated from 2 • 10 3 simulated MW-type galaxies. The contours are spaced 
by 1 M pc -2 for HI and 20 M Q pc -2 for H2, with the bottom contour cor- 
responding to the separation value. The redsh ifts z = 1.5 and z = 3 respec- 
tively correspond to the galaxy BzK-21000 (Daddi et al. 2008) and to the 
primary science goal of ALMA, i.e. the detection of MW-type galaxies at 
z = 3 in less than 24 hours observation. 



Redshift z 
Type 

M s /[10 10 M o ] 
M ffl /[10 9 M o ] 
rhf/[kpc] 
£hf/[M oP c- 2 

M H2 /[10 9 M o ] 
rj£ lf /[kpc] 

Ehalf/[M 0P C- 2 



Sim/MW 



Sim/BzK 



Sim 



0.0 

Sb-c/SBbc 



5+0.3 

7^-4/ 8 +2(b) 
16!^/15 + 5(b) 

A V6 +3(b) 
2 /3.5 + |( c ) 

2/7+1 (c) 
2 // -l 



1.5 

Sc/S? 
L7 + _° Q 2 2 /6 + _ 
2.3 



4.5; 
3.5; 

6.8 
12+ 



+07i 



7.3+ 
5.4 + 
4.6+ 
2.9+ 



7? 
/? 
/? 

/17+j 
/4!l 



130^0/150+50 



3.0 

Sc-d 

0.5+ 01 

7+^ 
3.9+ " 
5.5+ 

2.9+^ 

1 4+#3 

-0.2 
360^90 



Table 2. Average value and 1-cr scatter of the cold gas and stellar properties 
of 2- 10 3 simulated MW-type galaxies (Sim) at z = 0, z = 1.5, and z = 3. The 
actual measurements of the MW with 1-cr uncertainties were adopted from: 
(a)lFlynnetalJ d2QQ6h : (b) from anal ytic fits to Z ffl (r) inlKalberla & Dedesl 
(2008); (c) from S H2 (r) in Table 3 in landers etail dl984|) . The properties 
of the simulated MW-type galaxies at z =1.5 are c ompared to those of the 
CO-detected galaxy BzK-21000 daddi et al. 2Q0|)- 



5 CONCLUSION 

We have studied the cosmic evolution of the surface densities of 
HI and H 2 in regular galaxies using theoretical models combined 
with the Millennium Simulation. 

A key result is that the surface density of HI remains approx- 
imately constant and close to saturation at all redshifts, while the 
mean surface density of H 2 changes dramatically as oc (1 + z) 2A , 
mainly due to the size evolution of galaxies. These predictions will 
become testable with future telescopes, such as the LMT, ALMA, 
and the SKA. The few CO-detected high-z galaxies available today 
seem to have even higher H 2 - surface densities than predicted by 
the oc (1 + z) 2A scaling (see Fig. [2]). This could result from a selec- 
tion bias towards systems, which have lost some of their specific 
angular momentum in major mergers. 

We also studied the cosmic evolution of the cold gas in a 
sample of MW-type galaxies. We predicted that at z = 3 the 
H 2 -mass of a MW-type galaxy is M H2 ~ 3 • 10 9 M o with the 
denser 50% of this mass reaching an average surface density of 
500M o pc- 2 . 



Framework Programme 6, Square Kilometre Array Design Studies, 
contract no 011938. The Millennium Simulation databases were 
built as part of the activities of the German Astrophysical Virtual 
Observatory. 
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